Collective modes in uniaxial incommensurate-commensurate systems 

with the real order parameter 



V. Dananic 

Department of Physics, Faculty of Chemical Engineering and Technology, University of Zagreb 

Marulicev trg 19, 10000 Zagreb, Croatia 

A. Bjelis and M. Latkovic 
Department of Theoretical Physics, Faculty of Science, University of Zagreb 
Bijenicka 32, 10000 Zagreb, Croatia 

The basic Landau model for uniaxial systems of the II class is nonintegrable, and allows for 
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dimerized) ordering. In the present paper we complete the analysis of this model by performing 
the second order variational procedure, and formulating the combined Floquet-Bloch approach to 
the ensuing nonstandard linear eigenvalue problem. This approach enables an analytic derivation 
of some general conclusions on the stability of particular states, and on the nature of accompanied 
collective excitations. Furthermore, we calculate numerically the spectra of collective modes for 
all states participating in the phase diagram, and analyze critical properties of Goldstone modes 
at all second order and first order transitions between disordered, uniform and periodic states. In 
particular it is shown that the Goldstone mode softens as the underlying soliton lattice becomes 
more and more dilute. 

05.10.-a, 05.70.Fh, 64.70.Rh 



I. INTRODUCTION 

One of the most useful insights into the properties of stable and metastable ordered states in many body sys- 
tems follows from the investigations of accompanying collective modes, excitations with a coherent participation of 
(semi)macroscopic number of particles. The attention is usually focused on the lowest branch in the spectrum. If 
it is of Goldstone type, i. e. gapless (e. g. acoustic) in the long wavelength limit (k — ► 0), there is a continuous 
degeneracy in the characterization of ordered state, associated with the breaking of symmetry of the high temperature 
thermodynamic phase. Without the continuity in the degeneracy one has instead a finite gap at k = 0. 

Obvious extrinsic causes for the gap in the Goldstone mode are impurities, defects in the crystal structure, etc. 
Another cause for the gap is the presence of long range interactions jl| . Wc do not consider either of these mechanisms 
here, but remind that, as is well known in charge density wave materials [§J|, they may play a decisive role in the 
collective dynamics of ordered state. Instead, we concentrate on the systems in which the above distinction regarding 
the degeneracy of ordered state(s) has its origin in short range interactions. Those are numerous materials that show 
one or more types of uniaxially modulated orderings with periodicities which may be commensurate or incommensurate 
with respect to the underlying crystal lattice [f|J|] . 

Let us at the beginning invoke some simple widely accepted conclusions, accumulated through intense theoretical 
and experimental investigations on these incommensurate-commensurate (IC) systems in last few decades. In an 
ideal case of sinusoidal modulation the spectrum of collective excitations contains two types of modes, phasons and 
amplitudons, representing linearized fluctuations of phase and amplitude of the order parameter respectively. While 
the amplitudon mode has a finite gap below the critical temperature, the phason mode is acoustic if the free energy of 
corresponding state does not depend on the relative phase of ordered modulation and crystal lattice. In other words 
one has the continuous degeneracy with respect to this relative phase. It is strictly fulfilled only if the modulation is 
incommensurate with respect to the periodicity of crystal lattice. 

For commensurate modulations the free energy depends on the relative phase, as is easily seen already from the 
standard Landau expansions in which the lattice discreteness is taken into account by keeping a leading Umklapp 
contribution. Within this standard and frequently explored model ^,0], which leads to the simple variational equation 
of sine-Gordon type, the phason mode acquires a gap which is finite only for the strict commensurate ordering, and 
diminishes rapidly (exponentially) as the order of commcnsurability increases. For other modulations, which may 
have the form of dilute soliton lattices, the Goldstone mode remains gapless, although among these modulations there 
are solutions with commensurate periodicities close to the exempted leading commensurability. In other words, within 
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this model one does not distinguish between "secondary" commensurate orderings and incommensurate orderings. 
This is the consequence of a crude simplification made by retaining only one Umklapp term in the free energy. The 
recent analysis shows that already after taking into account two leading Umklapp terms the phase diagram becomes 
qualitatively different (|||. ^ contains a finite number of commensurate states, and shows a harmless staircase, i. e. 
a series of first order transitions between neighboring states. The Goldstone mode is then expected to have a finite 
gap for each state participating in the phase diagram. 

The Landau models for the orderings with spatial modulations are generally justified providing the interactions 
responsible for their stabilization are weak enough, so that the variations of order parameter (defined with respect to 
the appropriately chosen star of wave vectors) are slow at the scale of lattice constant. Two crucial simplifications 
are then allowed, namely the gradient expansion and the perturbative treatment of lattice discreteness through the 
truncation of the sum of Umklapp contributions. 

In the opposite regime of strong couplings the above spatial continuation is not allowed, and the lattice discreteness 
leads to qualitatively different properties of phase diagrams and related spectra of excitations, established by numerous 
analytical, and particularly numerical, studies of spin (e.g. Ising) Jf^,|ll[, displacive (e.g. Frenkel-Kontorova) [|||, and 
electron-phonon (e.g. Holstein) jl3],[l4| discrete models. Characteristically for such models, either a finite, sometimes 
large, number of commensurate modulations in the cases of harmless staircase, or the infinity of them in the cases of 
complete devil's staircases, can participate in the phase diagram. All commensurate states then have lowest branches 
of collective excitations with finite gaps in the limit k — ► 0. We repeat that, while none of these possibilities can 
be reproduced by Landau model with one Umklapp term, the former harmless staircases with a finite number of 
commensurate states are realized already within extended Landau models with only two Umklapp terms taken into 
account (|||. 

The analysis of Frenkel-Kontorova and Holstein models established also a new type of instability that involves 
incommensurate modulations, the so-called transition by breaking of analyticity JlSj . Namely, by increasing the 
coupling constant or by decreasing temperature jlTjJl^], the smooth envelope of an incommensurate periodic 
modulation becomes nonanalytic. The free energy then depends nonanalytically on the relative phase of modulation 
and underlying lattice. As a consequence a finite gap opens in the Goldstone branch of collective excitations even for 
incommensurate modulations. 

Already from the beginning of investigations on discrete models it was realized that the above complex features 
in phase diagrams and spectra of collective excitations have their origin in the nonintegrability of these models, i. 
e. in the nontrivial chaotic structures of corresponding phase spaces. In this respect it is important to emphasize 
that, cither in their basic form or after the inclusion of further terms, Landau free energy expansions are as a rule 
the examples of nonintegrable functionals. For example, while the sine-Gordon model with one Umklapp term, as 
a basic model for class I of IC systems, is integrable, already the inclusion of another Umklapp term brings in the 
nonintegrability ||PJ- 

The situation is even more intriguing for the class II, i. e. for IC materials with modulations having the period close 
or equal either to the original or to the dimerized unit cell of crystal lattice. There are numerical jl!| and analytical J2(| 
indications that already the minimal |2l|]22|l , as well as slightly extended HH , models for this class are not integrable. 
The consequences of this nonintegrability on the phase diagram are discussed in detail in Ref. po|| . In particular, it is 
shown that, in addition to simple disordered, commensurate [i. e. (anti)ferro] and (almost) sinusoidal incommensurate 
states, included into previous analyses |2^,^J,^3|, the phase diagram contains also an enumerable family of metastable 
solutions with the periodic alternations of commensurate and incommensurate sinusoidal domains. 

In the present work we calculate the spectrum of collective modes for stable and metastable states in systems of 
class II. The corresponding Landau model is particularly convenient for the discussion of questions raised in this 
Introduction, since it is nonintegrable, and, in addition, the accompanying phase diagram comprises both commen- 
surate and incommensurate (meta)stable states. Our main aim is to investigate to what extent are the collective 
modes influenced by the nonintegrability of, here continuous, free energy functional. Furthermore, by analyzing Gold- 
stone modes for the modulated states of the model under consideration |^2|,^3| we also resolve some controversies 
present in the literature p5|-p9t on its applicability in the description of incommensurate phases in systems of class 
II. The equivalent analysis for the class I, i. e. for the Landau model with two Umklapp terms, will be presented 
elsewhere p[| . 

The plan of the paper is as follows. The free energy functional for the class II is introduced in Sec. II. In Sec. Ill 
we perform the variational procedure up to the second order, taking care about some specific questions related to the 
thermodynamic minimization pf[ | . The linear eigenvalue problem associated to the second order variational procedure 
is discussed in Sec. IV. Here we encounter a generalized Hill problem, since the systems includes four coupled first 
order equations (in contrast to the standard cases with two equations), and furthermore, since we are looking for the 
collective modes of highly multiharmonic periodic states. We therefore do not follow a standard way, appropriate for 
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simple sinusoidal incommensurate orderings, but develop for the first time a general formalism, applicable also to other 
types of Landau models. This formalism enables the determination of Floquet exponents, and of corresponding Bloch 
basis of eigenfunctions which we consider in Sec.V. The numerical results for the collective modes of all (meta)stable 
states appearing in the phase diagram are presented in Sec. VI. Concluding remarks along the lines specified in the 
previous paragraph are given in Sec. VII. 

II. MODEL 

The free energy functional for the uniaxial incommensurate systems of class II is given by 



L 



dz (1) 



where u represents the real order parameter and L is the length of the system. This functional is the simplest 
(minimal) Landau expansion for the systems with minima of free energy density in the reciprocal space close to the 
center of Brillouin zone, or to the part of its border perpendicular to the uniaxial direction. Then c < 0, and one has 
to add a highly nontrivial term with the second derivative of u (and presumably positive coefficient d) in order to 
ensure the boundness of Landau expansion in the reciprocal space. The rest of the expansion (]]]) is standard, with 
b > 0, and a becoming negative below the critical temperature of the transition from disordered to uniform (ferro) 
or dimerized (antiferro) phase. We limit the further analysis to the most interesting regime characterized by c < 0. 
It includes the incommensurate ordering and the transition to the commensurate ordering (but does not include the 
transition from disordered to commensurate state which takes place for c > 0) pl| , ^| . In this regime the useful 
dimensionless quantities are 

L=J--L, u (z) = -^—u(z), f[u] = - r f[u}. (2) 
The model (0) can be now represented as the one-parameter problem, 



—L 



d 2 u\ 2 ( du\ 2 



dz 2 J \dz J 



) +\u 2 + ±u 4 



dz, (3) 



with A = ad/c 2 . The parameterization of the phase diagram in the regime c < is thus very simple, since all 
relationships between different (meta)stable states (like phase transitions, ranges of coexistence of two or more states, 
etc) can be presented in the one-dimensional A-space. The knowledge of the actual dependence of this parameter, as 
well as of the scales which enter into the reduced quantities (Q), on the original physical parameters, in particular on 
temperature, goes together with the specification of microscopic background behind the phcnomcnological free energy 
([!]) . This is a necessary step in any comparison of phase diagram for the model (|^) with experimental data for a given 
material. 

In the previous works Jl^,^] on the functional (||) we have determined thermodynamically stable states, i. e. 
its local minima, without taking into considerations statistical fluctuations outside these minima. This mean-field 
type of approximation is inappropriate for (quasi) one-dimensional systems. It is however usually sufficient for three- 
dimensional uniaxial systems with strong enough couplings in the perpendicular directions, on which we concentrate 
here. 

The thermodynamic extremalization of functional ([|) consists of the standard variational procedure that is equiv- 
alent to the classical mechanical one and leads to the corresponding Euler-Lagrange (EL) equation 

d 4 u d 2 u o 

^ + ^ +Au + U= °< (4) 

and of the extremalization that involves boundary conditions or some equivalent set of parameters. The general 
procedure that carefully takes into account the latter aspect is proposed in Ref. J3l|]. The most interesting result of 
this approach is obtained for the functional with the kernel that is not explicitly z-dependent. Then the relation 

F = -H, (5) 
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holds for each thermodynamic extremum uq(z). Here F is the corresponding averaged free energy, and H is the 
integral constant of the problem (JJ) which corresponds to the Hamiltonian in classical mechanics. 

Early considerations of the model (|l|) led to the suggestion that the mean-field phase diagram Sj^^jjEgLHI 
contains only disordered [ud(z) — 0], commensurate \u c (z) = ±y/—X], and (almost) sinusoidal [u s (z) as 
2/\/3(-v/l/4 — A) sin (z/y/2)] incommensurate orderings. The commensurate state is thermodynamically stable for 
A < —1/8 (and for a < in the range c > 0). The incommensurate state is stable in the range —2 < A < Aid = 1/4, 
while the first order phase transition between the commensurate and the incommensurate states occurs at Ai C = 
— 1.112. The more precise values, obtained after taking into account corrections from higher harmonics in the sinu- 
soidal ordering poj, are —1.835 < A < A^ and A; c = —1.177. Also, the wave number of this ordering, q, slightly 
deviates from l/v2 [i. e. from y—c/ (2d) in the original parameters of Eq. (jl])] as one approaches the left edge of 
instability, A — > —1.835. 

While by above solutions of EL equation (|4|) one exhausts all absolute minima of the free energy ([!]), the more 
involved numerical analysis [l9|^o[| showed the existence of an enumerable series of periodic solutions which arc 
metastable in finite ranges of the parameter A. The corresponding phase diagram is shown in Fig. [l] in which we 
ascribe to various solutions symbolic words introduced in Ref. ]|(| . By their physical content the metastable solutions 
from Fig. [j] represent periodic trains of successive sinusoidal and uniform segments (see Fig. 1 in pp| ), and, as domain 
patterns, complete in a natural way, as an inherent outcome of nonintegrable model (|l|), the phase diagram in the 
range of coexistence of two corresponding basic types of orderings. 



III. SECOND ORDER VARIATIONAL PROCEDURE 



The question on which we concentrate now is the thermodynamic stability of a given state u(z) which is a solution 
of EL equation ([I]) and fulfills additional conditions of Ref. ]3l| . To this end we have to go beyond the linear terms 
in the extremalization procedure. Let us therefore at first extend the standard variational procedure to the second 
order. Later on we shall shortly consider the conditions which follow from the minimization of boundary conditions. 

Let rj(z) be the infinitesimal variation with respect to u(z), obeying usual conditions at the boundaries z = and 
z = L, 

r)(z = 0) = rj{z = L) = rj'(z = 0) = rf(z = L) = 0. (6) 
After performing standard partial integrations, 

L J ^" >2dz = J 71 ' 71 ^ ~lJ 11 " lldz ' 

\ / V) 2 ^ = y'rit - ^'"'?lo + ~ [ V IV vdz, (7) 
the quadratic contribution to the corresponding variation of free energy functional (|^) can be expressed in the form 

= /[« + rj] _/[«] = I jf dz V V v , (8) 

with 

The linear differential operator (^) defines the eigenvalue problem 

V VA = V ^(z) + r}'i(z) + [A + 3u 2 (z)] VA (z) = A VA (z), (10) 

with the boundary conditions for t](z) specified by Eq. (JsJ) - The necessary condition for the thermodynamic stability 
of solution u(z) is given by the requirement that the spectrum A is non-negative for all normalizable solutions rj A (z) 
of the problem (^,||). 

Since the above procedure strictly respects the boundary conditions (|^), it is entirely equivalent to that usually used 
in classical mechanics. As a consequence the obtained condition for the stability of a given solution u(z) holds for any 
value of the sample length L. However, neither the extremal solution u(z) of EL equation (H), nor the corresponding 
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conditions of thermodynamic stability, should be sensitive to the conditions imposed on the sample surfaces in the 
physically relevant thermodynamic limit L — > oo. Therefore the stability condition can be generalized in this limit. 
In particular, we may ignore the boundary conditions (^|) , and perform the variational procedure for any infinitesimal 
variation rj(z), noting, for later purposes, that the requirement of infinitesimality excludes variations r](z) which would 
scale as |z|' 3 with (3 > 0. 

Performing the same steps as before, but now with neglected surface terms in Eqs. ^ (which scale as we 
come again to the linear eigenvalue problem (flO|), but without a specification on boundary conditions. This means 
that any complete set of eigenfunctions tja (z) with the eigenvalues A [which are themselves characterized solely by the 
linear differential equation flio|)] can be used in the representation of a given variation rj(z), and in the corresponding 
diagonal representation of the free energy (||). The condition A > thus guarantees the thermodynamic stability of 
the solution u[z) with respect to any infinitesimal variation, provided the system has the well-defined thermodynamic 
limit as specified above. We note that by this relaxation of boundary conditions @ we extend the standard "classical 
mechanical" second order variational procedure by including a part of, but still not all, "thermodynamic" variations. 
The discussion of this question in Appendix A suggests that the criterion of thermodynamical stability is probably 
entirely covered by the eigenvalue problem 

The concise definition of the thermodynamical stability, i. e. of the stability of any (absolutely stable or metastable) 
local minimum of thermodynamic functional with respect to small fluctuations, is thus: 

• A given solution of EL equation (ji|) is thermodynamically stable if and only if all solutions of Eq. ( |l0| ) for any 
A < are non-normalizable. 

For later purpose it is appropriate to introduce here also the concept of orbital stability, relevant for the behavior of 
particular solutions in the phase space: 

• A given solution u(z) of EL equation ([|) is orbitally stable if and only if all solutions of Eq.([lO|) for A = are 
normalizable. 

In the next section the above definitions will be used in the study of stability of homogeneous and periodic con- 
figurations u(z). The crucial assumption in this respect is that the solutions of Eq.([lO]) smoothly depend on both 
parameters A and A. 

Before embarking into the calculation of the spectrum of eigenvalue problem (|l^) , we invoke its general property 
which follows from the fact that the density of free energy functiona l (|3[ ) does not depend explicitly on the spatial 
coordinate z. Then there exists a normalizable solution of equation (noFwith A = 0, namely t)q{z) oc u'{z). This is 
the Goldstone mode that follows from the translational invariance of free energy functional ( |3]), by which u(z + zq) 
with arbitrary z is the solution of EL equation (Q) if u(z) is its solution. We note that Eq.([l0|) then has, together 
with the above Goldstone mode, another solution of the form rji(z) = w(z) + z ■ u'(z), where w(z) is some periodic 
function of the same period as that of the Goldstone mode. Although rji(z) is non-normalizable i.e. its norm grows 



as a power of L, we consider this non-normalizability as marginal. The power law growth of a solution of Eq.(10) 
is much easier to control than possible exponential growth of the remaining solutions, if there are any. For special 
values of the parameter A figuring in Eq.([l0|) with A = the only normalizable solution is the Goldstone mode u'(z), 
while other solutions have a power law growth in z, z n u'{z) with n < 3. These special values of A denote the edges 
of thermodynamical metastability of the corresponding configuration u(z). 



IV. FLOQUET THEORY 

The analysis of the eigenvalue problem ( |l~0| ) with periodic functions u(z) is based on general Floquet and Bloch 
theorems for linear differential equations with periodic coefficients. It will be performed in two stages, covered by this 
and the next section. The aim of the first one, based on the Floquet's approach, is to answer the question: Whether 
there exists a normalizable solution rj\(z) for a given value of A? In the second stage we calculate the set of values 
A for which normalizable solutions exist, i.e. the spectrum of collective modes, by using the Bloch's wave number 
representation. 

We begin by showing that the set of values of A for which the corresponding normalizable solutions T]\(z) may exist 



is bounded from below. To this end let us rewrite Eq. (10) in the form 



P J //\Ui^«Ur/n(:] ( A + ± - x\ VA (z) , V = -jL + I, (H) 
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and introduce the norm of the function t]a(z), 



lI 



|| VA || 2 = (tiXva) = y I VA(z)*VA(z)dz. (12) 



After multiplying Eq. ( |ll| ) by rj\{z) and integrating with respect to z we get 



|| V VA (z) || 2 +3 || u(z) VA (z) || 2 = (A + ~ - Aj || r? A (z) || 2 . (13) 

Here it is taken into account that the operator T> is hermitean and the function u(z) is real. Since the left-hand side 
of Eq. (|l|) is strictly positive, we conclude that 

A > Knin = X ~J ( 14 ) 

for each A for which the the norm ( |l2|) of the function T]a{z) exists. In particular, this means that it is sufficient 
to reduce a (numerical) analysis of the thermodynamic stability of a given configuration u(z) to the search for the 
normalizable eigenfunctions rj\(z) in the finite interval of A, A m i n < A < 0. 

Before considering Eq. jic| ) with the general periodic function u(z), let us establish the criterion for the thermody- 
namic stability of the particular homogeneous (jerro or antiferro) solution u c (z) = ±\/— A of EL equation (^J). Then 
Eq. ([to]) reduces to the linear differential equation with constant coefficients, so that the normalizable eigenfunctions 
must have the form f](z) oc e lkz with real values of the wave number k. The corresponding eigenvalues A are given by 

A = /c 4 — fc 2 — 2A , A < 0. (15) 

It follows that the homogeneous configuration u c (z) = ±V— A is stable, i.e. that A > for any k, provided that 
A < —1/8. Note that the latter inequality is just the condition of orbital instability of the homogeneous solution. 
Namely, the linearization of the EL equation with respect to this solution leads to the linear equation 

6"" + 6" - 2X0 = 0, (16) 

which has normalizable solutions 9(z) only for A > —1/8. Thus, we see that in this simple case the thermodynamic 
stability excludes the orbital stability, and vice versa, and that two stabilities "meet" each other in one point, A = — 1/8. 



A. General Floquet's procedure 



In order to apply the well-known Floquet's procedure |m| to Eq. ( |To| ) with a general periodic function u(z) 7 we 
rewrite this equation in the matrix form 



d&(z) 

dz 



A(z;X,A)&(z), 



(17) 



where ®(z) = [r)(z), r)'(z), r)"(z), r)'"(z)] , and the matrix A(z; A, A) is given by 



A(z) 



( 10 0^ 

10 

1 

y A - A - 3«(z) 2 -1 ) 



(18) 



The system of linear equations ( p"7j ) has four linearly independent solutions, &i(z),i = 1,...,4. They form the 
fundamental matrix 



F(z) = [® 1 (z),®2(z),® 3 (z),® A (z)}, 
which is obviously the solution of equation 

dF(z) 



dz 



A(z;A,A)F(z). 



(19) 



(20) 



G 



Without reducing generality we can always choose such initial conditions at z = that F(0) = I, where I is the 
identity matrix. 

The Floquet's theorem states that whenever the matrix A(z; A, A) is a periodic function of variable z with a period 
P, the fundamental matrix has the form 

F(z) = G(z)e Sz , (21) 

where G(z) is a matrix which varies periodically with z, G(z + P) = G(z), and S is a constant matrix. Due to the 
periodicity of matrix G(x) and the initial condition G(0) = I, the matrix S can be expressed in the form 

S = ilnF(P). (22) 

The matrix F(P) is called the monodromy matrix. The eigenvalues of £, c.j are Floquet's exponents and the 
eigenvalues of monodromy matrix F(P), p.j, are Floquet's multipliers. Put in other words, Floquet's theorem states 
that for each Floquet multiplier pi there exists a solution &i(z) of Eq. ( fi"7j ) with the property 

® z (z + P) =p&i(z). (23) 

Floquet multipliers pi are in general complex numbers. For a normalizable solution &i(z) the Floquet multiplier 
Pi lies on the unit circle in the complex p-plane, and the corresponding Floquet exponent Oi is imaginary. 



The problem (|17j) has an additional important property. 
[Z x {z), Z 2 (z), Z 3 (z), Z 4 (z)] T = T0(z), defined by 



B. Poincare-Lyapunov theorem 

After the linear transformation Z(z) 



Z x = i](z) , Z 2 = rj"(z) , Z 3 = 2(r/(z) + rj"'(z)) , Z 4 = 2rj'(z), 



(24) 



i. e. by 



Eq. ( p"^ ) is transformed into the equation 



T(z) 



dZ{z) 

dz 



/ 1 \ 

10 

2 2 

\ 2 J 



J H(z; A,A)Z(z), 



(25) 



(26) 



with 



H(z) 



/ -2(A- A + 3v (z) 2 ) \ 

2 

-| 

V -| i J 



/0 -1 \ 

-1 

10 

\0 1 / 



(27) 



The problem ( p6[ ) has the Hamiltonian form, characterized by the hermitean matrix H and the simplectic matrix J 
(i.e. J is antisymmetric and has the property J 2 = —I). The Poincare-Lyapunov (PL) theorem p5| ] for such problems 
states that the corresponding fundamental matrix, &{z), satisfies the relation 

$ T (z)J*(z) = $ T (0)J$(0). (28) 

In other words, $ T (z)J$(z) is the "integral of motion" for the Hamiltonian problem (|26|). This theorem can be easily 
checked by differentiating the relation (E8I) with respect to z, and taking into account Eq. (|2q). Since &(z) = TF(z), 



and the matrix T JT = Ji is also simplectic, it follows that 



F T (z)JiF(z) =F T (0)JiF(0) = Ji 



(29) 
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i.e. the PL theorem holds for our original fundamental matrix F(z) as well. 

From the relation (|||) it follows that the matrices F T (z) and F _1 (z) are similar [F T (z) = JiF _1 (z)J^ 1 ]. Thus, if 
p\ = p is the Floquet multiplier of F(P), then p^ 1 is also its Floquet multiplier. Furthermore, since in our example 
the matrix F(P) is real, it follows that p* and p*^ 1 are Floquet multipliers as well. These simple relations link four 
Floquet multipliers of the problem (|l7|) for any periodic solution of EL equation ([!]) and for any value of parameter 
A. The corresponding three possible types of distributions of Floquet multipliers in the complex p-plane are shown in 
Fig. |^. Floquet multipliers are either complex (a) or real. In the latter case two pairs generally have different values 
and may be of the same (b) or opposite (c) signs. Fig. || does not include the situations with existing collective modes, 
i. e. when one or two pairs of solutions are normalizable, and the corresponding Floquet multipliers are on the unit 
circle. 

Four Floquet multipliers of the problem ( p"7j ) can be represented as roots of a polynomial function of fourth order. 
Since, due to the PL theorem, p, p~ x , p* and p*~ x are all roots of such function, its general form is 

P 4 (p) - P 4 + o(A, A) (p 3 + p) + b(X, A)p 2 + 1, (30) 

where a(A, A) and 6(A, A) are, for a given periodic function u(z), some smooth real functions of parameters A and A 
from the matrix ([l8j). 



C. Scenarios of thermodynamic (in)stabilities 



Having recapitulated the Floquet theory for the Hamiltonian linear problem (17), we address the problem of 
thermodynamical stability for a given periodic configuration u(z). We start by noting that Eqs. (|l7j), ( |l8| ) and ( J30| ) 
enable some general conclusions about the dependence of the positions of Floquet multipliers in the complex plane 
on the parameters A and A. At first, since p — cannot be the root of Pi{p), it follows that by changing continuously 
A and A one can come from the distributions (a) or (b) to the distribution (c) in Fig. || only by passing through the 
unit circle. Next, it is easy to determine the positions of Floquet multipliers in the limits A — * — oo and A — > oo, 
since then we may neglect A and u 2 (z) with respect to A in the polynomial matrix element of the matrix A (|l|) [still 
keeping in mind that u(z) defines the period P which enters into the definition ((22" 

In the former limit A — > — oo the Floquet multipliers are given by 



fc " p , k n =| A | i e l< - 2n+1)2 i , n = 0,l,2,3 (31) 



Pn = C 

i.e. the distribution from Fig. |](a) is realized. Note that for A < A m j„ this distribution cannot pass to that from 
Fig. |^(c), since in this range of values of A the unit circle cannot be crossed because the problem ([l7]) does not have 
normalizable solutions. 

In the limit A — > oo the Floquet multipliers are given by 

p 1= p^=e^ p , p 3 = P 4 - 1 =e iAip . (32) 

As is seen in Fig. ||, one then has one pair of normalizable and one pair of non-normalizable solutions. 

From the other side, at A = one particular Floquet multiplier has the value p = 1, and corresponds to the already 
mentioned Goldstone mode. It has to be at least doubly degenerate, since otherwise the remaining three multipliers 
could not have symmetric positions required by the PL theorem. Possible distributions of Floquet multipliers for 
A = are shown in Fig. |]. Aside from the possibility that the degeneracy of the Goldstone mode is complete and all 
four multipliers are at p = 1 (a), one may have the remaining two multipliers either on the real axis (b,c), or on the 
unit circle (d). 

Taking into account the above conclusions, we are now able to list possible scenarios of thermodynamic (in)stabilitics 
for the periodic solutions of Eqs. (Q,|^). 

i) The solution u(z) is unstable for a given value of A if by increasing A from K m i n = A — 1/4 the Floquet multipliers 
from the distribution (a) or (b) of Fig. || move in such a way to come to the unit circle for some value of A in 
the interval (A - 1/4,0). 

ii) If the solution u(z) is thermodynamically stable, the Floquet multipliers for A — 1/4 < A < are defined either 
by one complex number not lying on the unit circle (Fig. |^a), or by two real numbers of the same sign (n, T2) 
and their reciprocals (Fig. |^b). The latter case has to be realized as A — > from below, since only distributions 
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(b) and (c) from Fig. [| represent the Goldstone mode for a thermodynamically stable configuration. Thus, at 
some negative value of A the distribution from Fig. ||a has to reduce to the double degenerate Floquet multiplier 
at the real axis (ri = 7-2), which then evolves into the distribution from Fig. ||b. 

iii) For special value(s) of control parameter (A = A c ) the thermodynamic instability of u(z) proceeds in a particular 
way, realized when all four complex Floquet multipliers approach together the point po — 1 as A tends to 
zero from below. The Goldstone mode is then completely degenerate (Fig. ^a). Putting in another way, such 
instability occurs when the points r and r _1 in Fig.^jb tend towards po = 1 as A — > A c . Note that the distribution 
of Floquet multipliers from Fig. ^c means that the instability, i.e. the crossing of Floquet multipliers with the 
unit circle, takes place at some negative value of A. Also, the distribution from Fig. [|d signifies that the 
remaining non-Goldstone mode is unstable in a finite interval of values of A, starting at some negative value of 
A. 



V. BLOCH THEORY 

In order to determine normalizable solutions of Eq. (|l7| ) with A > 0, i.e. the collective modes for given periodic 
configuration u(z) with the period P = 27r/Q, we profit from the freedom in choosing boundary conditions for the 
solutions rj{z), and specify periodic (Born - von Karman) ones. By this we chose the Bloch representation, 

% (z) = e* kz * k (z), (z+^) =* fc (z), (33) 



where k is the Bloch wave number limited to the I Brillouin zone (-Q/2 < k < Q/2). The differential equation for 
the periodic function ^k(z) reads 

^1 + Uk fl^l + (1 _ 6fc2 )^fl + 2 ik(l - 2^)^1 + 
az dz 4 dz dz 

+ [fc 4 - fc 2 + A + 3u(z) 2 ] V k (z) = A(fe)* fc (2), (34) 



and the normalizability condition is 



g- [ Q n{^k{z)dz = l. (35) 
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The dependence A(fc), i. e. the spectrum of eigenvalue problem (|f0|), follows from Eqs. |35| ). Since k is quasi- 
continuous in the limit L — > 00, this spectrum is for a stable configuration u(z) composed of non-negative bands. The 
corresponding Bloch functions "q n ,k{z), where n enumerates bands, represent a complete orthonormal set of functions 
for the problem (|l0|). 



From expressions (|23j) and (33) it follows that the Floquet multiplier for the Bloch function T]k(z) is given by 
p — e lhp . The polynomial function ( ^0| ) then has the form 

P 4 (p) = (p - e lkP ) (p - e- ifcP ) [p 2 + c 2 (A, A)p + l] , (36) 

where C2(A, A) is some coefficient. Comparing two representations for Pi{p) we conclude that the coefficients from 
the expressions (p0|) and (56|) are linked by relations 



a(A, A) = c 2 (A, A) — 2cos(/cP) , b(X, A) = 2 - 2c 2 (A, A) cos(fcP) , (37) 

i. e. that the eigenvalue A depends on the wave number k only through the function cos(fcP). This in particular 
means that for each band A(k) we have A(— fc) = A(k). This is consistent with the symmetry of Eq. (|34|). Furthermore 



A(k + = A(k), in accordance with the reduction of wave numbers in (33) to the I Brillouin zone. 



The representation (|33j) is particularly convenient for the analytical discussion of long wavelength limit k — > for 
the Goldstone mode for which A(k = 0) = and , 5a; = o(z) = u'(z). To this end we insert the Taylor expansions for 
small k, 

= u'[z) + + fc 2 * 2 (z) + ... , A(k) = k 2 A 2 + k 4 A A + (38) 
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into Eq.( |34 
equations 



and 



The requirement that the coefficients in front of leading powers, k and fc 2 , vanish then leads to the 



(39) 



(40) 



with the operator V given by Eq. (|TT|). After multiplying equation ( |40| ) by v', integrating with respect to z, using the 
fact that vl is the Goldstone mode, and inserting 2u"" + u" from Eq. (|39"|), we get the expression for the coefficient 
A 2 , 



V 2 + A - - + 3u 2 ) *! = -2i (2u"" + u") 



V 2 + A - - + 3m 2 ) * 2 = (A a + l)u' + 6u"' - 2i(2V'{' + 



A 2 = - 



1 + 6 7v2\ J7jT\ • 



(41) 



where (...) stands for the spatial integration, like in Eq. (|T^). However, the general thermodynamic condition (|J) for 
the functional (m reads [Rll 



(u" 2 ) 1 
(u> 2 ) ~ 2 ' 

so that Eq. (H|) can be written in a more transparent way, 

A 2 = 2 - F 2 [u{z)} . 

Here we introduce the functional 

(v 2 + A - | + 3u 2 ) *i 



^2 [«(«)] 



(42) 



(43) 



(44) 



Since the operator figuring in this equation just defines the eigenvalue problem ( |lO| , |ll[) , it is clear that the functional 
i*2 [u(z)\ is positive definite for any thermodynamically stable configuration u(z). This has two consequences. 

Firstly, the common upper limit of the velocity of Goldstone mode, vq — VA2, for all thermodynamically stable 
periodic states is vg,m = v2. The velocity vg for a given (meta) stable state has the maximum value va.max < vg,m 
when the functional (^ZJ) attains its minimum. Like the functional (^), the functional ^[^(z)] depends only on the 
parameter A. Thus, taken a given solution u(z), we can find the function VPi by solving the inhomogeneous linear 
differential equation (|39|), and then determine, by calculating F2[u(z)], the velocity vq as a function of A. In other 
words, we have a direct method for the calculation of the velocity of Goldstone mode, not related to the above 
Floquet-Bloch procedure (but derived from it). It can be used as an independent check of numerical results for the 
spectrum A(fc) which follow from Eq. (^i|). 

Secondly, the functional (44) attains its minimal value (F2 = 0) if and only i/thc function ^i(z) vanishes. As is 
seen from Eq. (|39|), this is possible only when the solution u(z) satisfies the equation 



1 "" 1 " n 

2u + u = 0. 



(45) 



i.e. when u(z) cx sin (z / y/2) . The only solution from the phase diagram in Fig. [I] with this property is the almost 
sinusoidal incommensurate state, denoted by s 2 . Since in the limit A — > Ay = 1/4 it reduces strictly to the above 
simple sinusoidal dependence on z (with the amplitude tending to zero), we conclude that just at the second order 
phase transition from the incommensurate to the disordered state the velocity of Goldstone mode of incommensurate 
state attains the maximum value vq.m = V%- We note that other periodic (and metastable) states u(z) from Fig. |l| 
cannot even approximately satisfy Eq . ([45|). From the other side, due to the deviations from the sinusoidal form 
of a given solution, the functional ( |44[) can attain the value F2 = 2, in which case the velocity of Goldstone mode 
vanishes. As will be seen from the numerical results in the next section, this is indeed the case for all periodic solutions 
(including the almost sinusoidal configuration s 2 ) at the edges of their local thermodynamic stabilities. 
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Let us conclude this general discussion with the remark on the class of periodic solutions u(z) which in addition 
have the property u(z + P/2) = —u(z). Since only u 2 (z), which then has the period P/2 (and not P), enters into 
the problem (|]j]), the corresponding spectrum A(fc) and the Floquet multipliers can be calculated with respect to the 
former period. The I Brillouin zone is then doubled (— Q < k < Q), and the number of branches of collective modes 
is halved. In particular, within this choice the value of Floquet multiplier for the Goldstone mode u'(z), defined by 
Eq. (p3|), is —1 and not 1. The approaching of A = from below for the stable solution u{z) then proceeds like in 
the point (it) of Subsection IV. A, but with one pair of Floquet multipliers tending towards the point po = — 1, and 
the other pair placed at the negative real semiaxis (Fig. ^c). Correspondingly, the Bloch representation of Goldstone 
solution of Eq. ( JTt| ) is r) = e ±l( 2 2 \l/(z) with = e Tl ® z u' (z), i.e. the Goldstone mode is placed at the border of the 
doubled Brillouin zone. However, the propagation of collective modes takes place in the periodic structure determined 
by the configuration u(z) (and the period P). Thus the above doubled Brillouin zone has to be folded once to get the 
physical one, —Q/2 < k < Q/2. In other words, the wave numbers k = Q and k = coincide, so that the Goldstone 
mode is realized as the long-wavelength one for such solutions as well. Furthermore, we note that after this folding 
the above Taylor expansion (|3^) and subsequent conclusions on the velocity of Goldstone mode [Eqs. 
in the same way for states with the property u(z + P/2) — —u(z) as well. 



VI. COLLECTIVE MODES FOR SYSTEMS OF CLASS II 

In order to derive collective modes for configurations participating in the phase diagram from Fig. lL we extend the 
numerical method developed in Refs. (lj]]2(]] to the calculation of eigenvalues and Bloch solutions (33) of the linear 
problem ( [To| , |l7| ). In the further discussion we shall mostly use the notation Q(k) = A(fc), where il(fc) has the 
meaning of the frequency of collective mode. Note that the energy scale for fi(fc) [as well as that for averaged free 
energies in Fig. |jj is defined by the last expression in Eq. (||). 

The periodic solutions from the phase diagram were determined by solving a system of algebraic equations for 
coefficients of their Fourier expansions. These Fourier sums were truncated at finite degrees, high enough to ensure 
a sufficient precision for u(z), as well as for corresponding wave number Q, averaged free energy, etc. The limitation 
of this method comes from the increase of number of non-negligible Fourier coefficients as the period 2n/Q increases, 
and the corresponding solutions u(z) contain more and more elementary sinusoidal and uniform segments. 

Representing the function u 2 (z) in Eq. (|^) by corresponding truncated Fourier series, and writing the function 
^k(z) in the same manner, 

N 

^k(z) = a + V2 2J [a n cos(nQz) + b n sin(nQ.z)] , (46) 
n=i 

we come to the homogeneous linear algebraic system for the coefficients do, di, .., fl/v and &i, &2> •••> &zv- I n order to 
calculate collective modes A(fc), it remains to diagonalize the corresponding (2N + 1) - dimensional matrix. This 
matrix is generally complex and hermitean. Again, one has to keep a sufficient number of Fourier components in the 
expansion (^6|) to get a reliable result for at least two lowest branches in the spectrum A(k). In actual calculations 
the truncation at a given number of coefficients N is taken as acceptable if for a given branch A(fc) one fulfills to a 
certain degree of approximation the equality A(fe = 0) = A(k = Q) [i.e. the equality A(k = 0) = A(k — 2Q) for the 
solutions with the property u(z + P/2) = — u(z)]. 



A. Collective modes of states u s (z), u c (z) and Ud(z) 

In Fig. ^|we present the spectrum of collective modes for the incommensurate almost sinusoidal state u s (z), denoted 
by s 2 in Fig. ^, choosing few characteristic values of the parameter A. As announced above, we use here the reduced 
Brillouin zone, —Q/2 < k < Q/2, for all values of A, except for those for which the periodic modulation is absent 
(A = 0.3 > \ ld in Fig. |c). At the very second order transition from the incommensurate state to the disordered state 
u d( z ) — (A = Ajd = 1/4) we present the spectrum in both, reduced and extended, zone schemes (Fig. |^c). Note 
that due to the additional symmetry of s 2 state, u s (z + tt/Q) — —u s (z), the subsequent branches in Figs. |^a-b are 
not separated by gaps at the zone edges k = ±Q/2. 

At first, we see that the lowest branch has the property of Goldstone mode [fi(fe) ~ k for k — > 0] in the whole range 
of stability of the configuration s 2 . For A well below the critical value Aid (A = —0.1 in Fig. ||a) the subsequent pairs of 
branches defined in such way are separated by gaps at k — 0. In other words, the general property obtained before in 
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the limit A = il 2 — > oo by which only one pair of Floquet multipliers is on the unit circle (Fig. ||) , is here realized for 
all values of A. However, as A increases the gap between two lowest pairs of branches decreases, and finally disappears 
for A w 0.05, as is seen in Fig. |a. Then one has an overlap of branches in a finite range of values of f2, i.e. all four 
Floquet multipliers are on the unit circle. This overlap increases, and the minimum of higher branch tends towards 
0, as A approaches the critical value Ay (A = 0.249 in Fig. |^b). At A = Ay this minimum has the value f2 = 0, while 
the slope d£l(k)/dk has a finite value which coincides with that of already existing Goldstone branch (Fig. ||c). In 
other words, just at the second order phase transitions one has two acoustic modes, which, although with same phase 
velocities v = lim^o > have different dispersions at finite values of k. For A > Ay these two branches combine 
into a single mode which has minima at k = ±Q with a finite value £l(Q), and a maximum at k = (Fig. |c), as it 
follows directly from the quadratic part of Landau expansion (|^). 

The dependence of the phase velocity of Goldstone mode, vq, on the parameter A is shown in Fig. [| It is finite 
at A = Ay, decreases as the amplitude of the incommensurate state increases, and vanishes at the metastability edge 
for the s 2 state, A = —1.835. This dependence is in accordance with the analytic results (41-45]) on the asymptotic 
behavior of the Goldstone mode. 

The spectrum of collective modes for the commensurate state u c {z) = [i.e. u c (z) = ±^/—a/b in the original 

notation of Eq. (|l|)] follows from Eq. (|l5|). This state is thermodynamically stable in the range a < for c > and 
a < —c 2 /(8d) for c < 0, which comprises positive values of c, excluded from the analysis after the transformation (0). 
In order to cover the whole range of stability of u c (z), we rewrite Eq. ((l^) in the original notation, 



O 2 = d~k 4 + c~k 2 - 2a = d ( k 2 + — ) - — - 2a (47) 

V 2d) 4d 



[with k = ^J—c/dk, Vl = ^/ 'c 2 JdO, in the range c < 0]. The second equality in the expression (|47]) shows that 
for c < the dispersion curve has minima at k = ±y / —c/(2d) (i.e. at k = ±l/\/2), with the gap tl(k) equal to 
y— 2a — c 2 / (4d) [i.e. y/2(— A — 1/8) in the reduced scale fl]. As for the range c > 0, it follows from the first equality 



in Eq. (15) that the collective mode has the minimum at k = 0, with the gap O(0) = y— 2a. The gap vanishes at 
a = 0, i.e. at the line of second order transition from the commensurate state u c {z) to the disordered state iid(z) = 0. 

The commensurate solutions u c (z) = ±y- a/b, which here represent the uniform or dimerized ordering for the 
Landau expansions ([!]) around the center or the border of the original Brillouin zone respectively, have the same 
symmetry properties as the solution for the disordered state, Ud(z) = 0. The only collective excitations with finite 
activated frequencies are fluctuations of the amplitude u{z) with the above dispersion relation (^). Note that the 
mode of Goldstone (acoustic) type is absent. Since the solutions u c (z) possess, as constants, a trivial translational 
degeneracy, we prefer to associate this absence of acoustic branch with its reduction to the trivial dependence Ct(k) = 0. 
The purpose of this interpretation will become clear in the next subsection. 

Finally, as it follows directly from the expression (Q), the disordered state Ud{z) = which is stable in the range 
a > 0, c > — V4ad, has a branch of collective excitations with the minimum at k = for c > 0, and with two minima 
at k = ±y— c/2d for c < 0. The respective gaps at these minima are equal to ^fa (for c > 0), and to \J a — c 2 /Ad (for 
c<0). 



B. Collective modes of periodic metastable states 

An illustration of spectra of collective modes for metastable states is shown in Fig. [j]. We take the state sd, chose 
the value of control parameter somewhere in the middle of corresponding region of stability from Fig. | (A = -1), 
and plot four lowest branches of collective modes. Spectra for all other metastable states from Fig. [j] have the same 
qualitative properties, and therefore are not plotted. More specifically, for all states, and for all values of A within the 
respective ranges of stability, the subsequent branches are separated by finite gaps, i. e. there is no branch overlap, 
like that obtained for the state s 2 [Figs. ||a,b]. 

Furthermore, the lowest branch for all states is the Goldstone mode with the dispersion fi(fc) ps vck for k — * 0, 
and with the corresponding phase velocity vg vanishing for the values of parameter A at the edges of stability. The 
dependence of vq on A for all metastable states from Fig. |l| is shown in Fig. M. The characteristic scales for these 
velocities, given by maxima vc.max of curves vq(X) for each metastable state, are situated in the range of values 
(0.4 — 0.5 in dimensionless units of Fig. |8|). This is to be compared with the maximum value of about 1.4 of vq for 
the configuration u s (z) (Fig. |^). In this respect one may recognize a rough tendency by which VQ^ max decreases as 
the proportion of incommensurate (s) domains decreases. This decrease is particularly evident as one compares s 
with s 9 d 3 , and with other configurations from Fig. 1 of Ref. EG], in which the proportion of commensurate domains 
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d is larger and both types of domains become rather short. The decrease of vo,max is then saturated, i. e. values of 
VG,max are roughly concentrated in the narrow range 0.47 — 0.48. 

The above tendency can be plausibly interpreted along the lines from the preceding subsection. The metastable 
states are in fact domain trains, built as successions of segments with local sinusoidal (u s ) and commensurate (u c ) 
orderings. As was already stated, it is plausible to associate to a commensurate segment a Goldstone mode with the 
vanishing frequency (and the vanishing velocity as well). The total Goldstone mode, which is some hybrid of these 
vanishing contributions and the contributions from the local sinusoidal ordering, tends to be softer and softer as the 
train has more and more commensurate domains. As a consequence, the velocity vc,max gradually decreases as the 
proportion of commensurate segments in metastable states increases. 



VII. CONCLUSIONS 



The results presented in Sec. 5 show that the spectra of collective excitations for all periodic states, stable and 
metastable, from the phase diagram of the model (@,|) [Fig. § have Goldstone branches with a linear dispersion 
fi = vck in the long wavelength limit. Thus, although these spectra belong to the nonintegrable model, they have 
standard characteristics that essentially follow from the absence of an explicit x-dependence of free energy density 
in Eq. (|l|). The latter property of the free energy in turn ensures the translational degeneracy of all solutions of EL 
equation (jij), including those participating in the phase diagram. In this respect the present spectrum does not differ 
qualitatively from those of integrable models with the same property. 

The fact that the chaotic content of the phase space for nonintegrable models, like for that defined by Eq. ([!]) Jl9| , p0| 
or for other examples |3^^|| , does not have as substantial impact on the spectrum of collective excitations as it has 
on the thermodynamic phase diagram, can be interpreted in the following way. The states from the phase diagram 
belong to the subset of solutions of EL equation defined by conditions like Eq. (||) . They are localized in the orbitally 
unstable chaotic layers which cover the phase space, have the measure zero in this space, and are mutually separated 



by topological barriers with characteristic heights given by averaged free energies of these layers |$7 38 1. These 
barriers do not allow for smooth changes from one state to others, and as such represent an intrinsic mechanism for 
frequently observed phenomena like memory effects and thermal hysteresis, as discussed in detail in Ref. |20[ . In 
general, nonintegrable free energy functionals have more complex phase diagrams than integrable ones. 

On the other hand, collective modes belong to another space of states, denser than the phase space, i. e. to that 
defined by the second order variational procedure and the corresponding eigenvalue problem (|l^). All states in this 
space are realizable as thermodynamic fluctuations. They have usual properties of double periodic linear systems, 
although the corresponding Bloch functions ^(z) in Eq. ( p3| ) may be far from a simple sinusoidal form. These 
properties are not essentially dependent on the level of integrability of free energy functional. 

In order to resolve the eigenvalue problem ( |Io| ) for the model (Q, ||), we formulate here a method based on the 
general Floquet-Bloch formalism, applicable to any IC system showing stable multiharmonic (i.e. non-sinusoidal) 
periodic ordering(s). Beside being a basis for the numerical calculations of eigenvalues and eigenfunctions (Sec. 5), 
this approach clearly indicates that for more complex models and orderings the traditional notions of phasons and 
amplitudons are not appropriate. In particular, it was often claimed that, being an expansion in terms of real order 
parameter, the functional (|l|) itself is insufficient for the stabilization of modulated states in the systems of II class, 
since incommensurate states, in particular those with soliton lattice like modulations, should have to be described by 
at least a two-dimensional order parameter |^,^,^,^9|. Also, the absence of phase variable in Eq.(|]) caused a belief 
p7| that the states which emerge from this functional do not have an acoustic (phason-like) collective mode. However, 
while the previous study [po| led to the conclusion that almost sinusoidal and highly non-sinusoidal configurations are 
among (meta)stable states of the model (0) (as is seen in Fig.p, the present analysis shows that Goldstone modes are 
well defined for all these configurations. From the other side, all dispersive modes for the homogeneous (u =const.) 
states are massive, i.e. have finite gaps. 

The gap of the lowest such mode tends to zero at continuous (2nd order) phase transitions from one homogeneous 
state to another, or to some periodic ordering. The examples are the lines (c > 0, a — 0) and (c < 0, A = ad/c 2 = 1/4), 
representing the transitions from the disordered state to the commensurate and incommensurate states respectively. 
As is shown in Fig. ||, the situation is qualitatively different at the transition from the disordered state to the 
incommensurate, almost sinusoidal, one. The reason is the specific behavior of Goldstone mode in the incommensurate 
state. By approaching the transition from the incommensurate side the phase velocity of this mode, vq, remains finite, 
while, as is shown elsewhere |59|, its oscillatory strength tends to zero. In fact, the above behavior of Goldstone mode 
for the state s 2 at the second order transition to the disordered state is exceptional. Namely, the Goldstone modes in 
the (meta) stable states behave critically at the edges of stabilities for these states, including the lower edge of state 
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s 2 at A = —1.835. At these edges the phase velocities vg vanish. All these specific properties of collective modes, 
particularly of the most interesting Goldstone modes, are expected to be directly experimentally observable in X-ray 
and neutron scatterings, as well as in optical and similar measurements. The particular discussion on the role of these 
collective modes in the dielectric response, and the comparison with measurements on some materials of the II class, 
is given in Ref. p9| |. 

Finally, we comment on the general property of Goldstone modes for metastable periodic states by which they 
become softer and softer as the period of these states increases. This tendency, shown in Fig. [§, has its origin in the 
elastic nature of Goldstone modes in the long wavelength limit. More specifically, as the segments of local sinusoidal 
order become more and more dilute in the underlying commensurate background, the slight variations in their mutual 
distances cost less and less energy, i. e. the corresponding effective elastic constant decreases. In this interpretation, 
which holds for dilute soliton lattices as well, the commensurate ordering is by assumption perfectly elastic, i. e. the 
notion of relative distance has no sense since the lattice discreteness is neglected. The only possible deformations are 
those invoking the variations of amplitude, and resulting in the massive collective modes. The lattice discreteness 
introduces, through an "external" potential of Peierls-Nabarro type, the finite stiffness of the local commensurate 
ordering, or even opens the gap in the Goldstone mode for dilute incommensurate states at the transition by breaking 
of analyticity 0011 • 
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APPENDIX A: 

The procedure from Sec. 3 takes into account, after relaxing boundary conditions (^), all infinitesimal variations 
of the order parameter u(z). This generalization includes some thermodynamic variations, like e. g. those, specified 
by the scaling u(z) — > su(z) with s — ► 1, responsible for the condition obeyed by u(z) at boundaries z — and 
z = L (condition B in Ref. pl|). However, by this procedure the analysis of thermodynamic stability is still not 
completed, since there remain variations which invoke infinitesimal relative changes in the configuration u(z), but are 
not infinitesimal at the absolute scale. An example is the scaling 

u{z) — >u[{l + e)z], e->0, (Al) 

which leads to the condition ([5]). The variation that corresponds to this scaling is not infinitesimal. Indeed, after 
the transformation (1 + e)z — > z in the integral (||), it follows that this variation behaves as z and therefore does not 
fulfill the criterion of infinitesimality specified in Sec. 3. Thus the above procedure has to be enlarged by including 
the expansion of the free energy with respect to e up to the quadratic terms. While the requirement that the linear 
term vanishes gives the condition (0), the second order variation reads 



S 2 f = f [«((1 + e)z) + r,] - f[u(z)] = 

' dz [r){z)Vr](z) + 2{u'(z)) 2 e 2 + A(2u""(z) + u" {z))r)(z)e\ , (A2) 



1 >' 
L 



i.e. the expression (g) is extended by the term quadratic in e, and the term representing the bilinear coupling between 
rj{z) and e. 



The previous analysis [g_9 20 1 of the model (|l]|3|) led to the conclusion that all solutions u(z) of the EL equation (^) 
that participate in the thermodynamic phase diagram as stable or metastable configurations, are simple periodic. The 
analysis in Sec. 4 shows that the corresponding eigenfunctions t]a(z) of the problem ( |ic| ) are then double periodic. 
This means that for periodic extrema u(z) the bilinear coupling in the expression ( |A2| ) vanishes, i. e. the fluctuations 
in e are decoupled from rj(z) fluctuations. The remaining e 2 term is positively definite, i. e. all periodic configurations 
satisfying the EL equation ([IJ) and the condition (||) are also stable with respect to the variation defined by the 
scaling (Al), irrespectively to the value of the control parameter A figuring in the functional (||). 
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FIG. 1. The phase diagram for the model (|3|). The inset shows the whole range of (meta)stability of the almost sinusoidal 
configuration s 2 and indicates the portion of the phase diagram in which other metastable non-sinusoidal configurations coexist. 
This portion is shown in the larger part of the figure. 
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FIG. 3. The distribution of Floquet multipliers corresponding to the existence of just one complex normalizable solution (a 
collective mode) of Eq.(|rij). 




FIG. 4. The completely degenerate Goldstone mode (a) which originates by merging two real multipliers shown in (b). The 
situation shown in (c) corresponds to approaching the instability for some negative A, which then evolves by moving two 
multipliers (corresponding to a non-Goldstone mode) along the unit circle. 




FIG. 5. The dispersion curves for the almost sinusoidal configuration for different values of the control parameter A. The 
curve (a) shows the situation when two branches are separated (A < 0.05). By increasing A they begin to overlap for A ~ 0.05 
(b), until A reaches the critical value A = 0.25. The curve (c) corresponds to A = 0.249. The finite slope of curves at k = 
coincides with that of the Goldstone mode. For A > 0.25 (d), these two brances combine into a single mode. 
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FIG. 7. Four lowest branches of collective modes for the simplest of non-sinusoidal state sd for A = — 1. The branches for all 
other non-sinusoidal configurations are qualitatively same as that for sd. 
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